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syhopsi s 


The present work is based on two separately developed 
concepts of procurement and inventory operations. For these 
concepts, two models are proposed. The first model deals 
with a multi-item inventory system involving single source 
and limited warehouse capacity. A dynamic programming 
formulation with usual cost terms is structured. Objective 
is to minimize the total cost per period. The total cost 
is expressed as a function of maximum inventory level. 

Depending on the situation under consideration, this variable 
can take discrete or continuous values. Therefore, both 
discrete and continuous dynamic programming have been used 
to solve the problem. 

The second model deals with joint replenishment of 
several items. The problem has '■ been formulated mathematically 
with the objective of minimizing the total cost. Three 
heuristic approaches have been developed. Heuristic one 
considers a planning horizon of finite duration. Heuristics 
two and three solve the same problem but the planning horizon 
is taken as infinite. For an n-item problem heuristic two 
generates n feasible solutions, and the feasible solution 
giving the minimum objective function value isselected. 
Heuristic three is an iterative procedure and it generates 
only one solution. 



viii 


The three heuristics have been compared for a set of 
45 problems generated randomly. It is found that heuristic 
three which is computationally very efficient yields better 
solutions in most of the cases. 



CHAPTER I 


INTRODUCTION 


The control and maintenance of inventories of physical 
goods is a problem common to all enterprises in any sector of 
a given economy. Inventories are maintained by manufacturing 
as well as non-manufacturing organizations. The economic 
significance of effective control of inventory is without 
question very great, for the costs associated with procure- 
ment and holding of inventory are large in almost all industrial, 
business and service organizations. It is no secret that 
inventories consume most of the investment capital in any 
enterprise. In fact inventories can represent more than 
50 percent of a company's invested capital and as much as 
80 percent of its cost of goods sold [lo]» Because of the 
magnitude of costs involved and the widespread applicability 
of inventory problems, it is obvious that inventory theory is 
a fertile field for cost reduction efforts. 

Al though literature on inventory can be traced back to 
the turn of the century, useful work has ensued only after 
1954. The later work attempts at structuring real life 
problems. The formulations, however, have been so complex 
and cumbersome that most of the research workers have concen- 
trated on only single product inventory models. 
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The present -work can he seen as consisting of two 
parts. The first part presents a model for developing the 
optimal procurement and inventory policies with respect to 
a cost function in a deterministic system composed of many 
items that may he procured from a specified source. The 
system is subject to a restriction on total warehouse 
space. The policy variables for the system under considera- 
tion include the procurement level and economic order 
quantity for each item. These are determined with respect 
to the relevant costs, system parameters, and a warehouse 
restriction such that the minimum total cost results. The 
minimization is performed through the application of dynamic 
programming technique. The solution proceeds stagewise 
(where each stage corresponds to an item) with the aid of 
recurrence relations and a functional equation method. 
Computer programs for both dynamic programming (discrete) and 
dyn am ic programming (continuous) have been developed, sasd 

* % 

The second part of the presentation aims at finding 
the optimal inventory policies (i.e. when to order and how 
much to order) with respect to a cost function when it is 
possible to replenish many items jointly. In case of joint 
replenishment the economic lot size for each item cannot be 
determined without considering the rest of the items within 
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the group. Bach time a hatch of product is procured,- there, 
is a common set-up cost for procurement in addition to the 
individual set-up cost for each of the items. Therefore, 
the economic order quantity for each item cannot he deter- 
mined without considering the rest of theitems which are 
jointly replenished with the item in question. 

The economic incentive behind thedesire for joint 
replenishment as a group lies in the fact that if each item 
is replenished individually, it must he responsible for the 
full set-up cost for each of its procurement run. In a 
manufacturing context, economies may accrue through joint 
replenishment as a consequence of savings in facility set-up 
and changeover time. In jointly ordering a number of items 
from a vendor econ omiesmay also accrue as a consequence of 
the preparation, processing and expediting of fewer orders 
and the receiving of fewer shipments. 

To solve such problems one common approach is to 
determine independently the lot size for each item as if no 
inter-relationship among items exists, and then to coordinate 
the replenishment of the items by modifying the lot sizes and 
cycle times of the items. One approach used in practice is 
to decrease the cycle time (and, accordingly, the order 
quantity) of all the items with a low order frequency, making 
it equal to the highest frequency item. The shortcoming of 
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this approach, however, is that it is unreliable. The 
method does not always yield a solution with minimum total 
cost. In fact, it sometimes yields solutions which arc 
more costly than those in which all cycle times are simply 
left at their original values. 

In the following' chapters, three different algorithms 
are presented for determining minimum cost solutions to 
multi item problems of this type. The first algorithm 
assumes a planning horizon of a finite duration and the 
last two algorithms assume a planning horizon of infinite 
duration. Solutions are sought in which all demands are 
met at exactly and forwhich a reasonably good value of objec- 
tive function is obtained. 



CHAPTER II 


LITERATURE SUR7EY 


Inventory control is a field which is rapidly growing 
in terras of available literature and am attempt to assess 
the work which, has been done already in this field and 
placing it in its proper perspective appears to be an uphill 
task. Keeping in view the objective of the present work, thi 
chapter concentrates on ’Multi-product Inventory Models'. 

An attempt has beon made to give an account of the state of 
literature as it stands today. 

Starr and Miller [23 ] appear to bo the first to 
introduce the concept of constrained multi-producb inventory 
systems. They assumed instantaneous replenishment and 
backorders wore not permitted. With the use of Lagrangian 
multipliers, an optimal stocking policy was determined for 
a given set of physical and economic constraints . Problem of 
this type has also been investigated by Hansman [ 7 ], Buchan 
and Koenigsberg [ 5 ], and Hadley and Whitin [ 6 ]. 

Haddor [L 7 ] considered a two product system and found 
two optimal policies. The first policy pertained to the 
joint ordering of products using the classical lot-sizo 
formula. In the second policy independent ordering of items 
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was permitted. Naddor and Saltzmann [18 ] considered a 
multiproduct system where the objective function was e sum 
of carrying and ordering costs. However* the ordering cost 
was a function of items listed on a single order. The prob- 
lem was solved using numerical approxim ations ana it gave 
subop timal results. Davis [ 4 ] modified the Maddor and 
Saltzman* s approach. He formulated an objective function 
that consisted of ordering* interest, obsolescence and 
holding costs. The holding cost tern included coots for 
warehouse spare, labour and equipment. However, since the 
obsolescence, interest, and holding costs wore all of the 
same form (a percent of the value of the item) they were 
consolidated into a single term that wo . 3 similar to carrying 
costs considered in previous in we s tiga tions . As before in 
the Had dor and Salzman model, the terms wore independent and 
the total cost was simply the sum of costs for each item 
considered independently. 

Plossl [20] proposed a LIMIT system to find optimal 
ordering policies, onco tho number of orders and reorder 
quantity aro known. LIMIT is an acronym for Lot -size Inventory 
Management by Interpolation Technique. It takes into account 
tho practical limitation on the total number of sot-up hours 
which are available. Optimal order quantity is calculated by 
the Wilscn's Lot Size formula and then the number of orders to 



satisfy the demand are determined. 3y an iterative proce- 
dure the optimal mode of ordering is found for a given multi- 
product inventory problem. The method was designed to reduce 
the total ordering costs. The method, however, did not 
recognise smy product interactions in the form of limited 
storage space or budgetory restriction. 

Shu [22] has described a method for determining 
replenishment frequencies of two items jointly replenished. 
Essentially, the procedure developed a shipping rule for 
items. The skipped item is replenished every E-th time (E 5^1) 
that tho orders are placed. A set of graphs is provided to 
determine E and H (where N is the number of times tho product 
will be procured). If the number of items exceeds two, Shu 
suggested that the items be ranked according to decreasing 
order of sales and that the smallest demand item be chosen 
as skipped item. The same procedure is to be repeated until 
the value of E for the item in question happens to be 1 or 
there are no more items to evaluate except the largest demand 
item. Earlo [ 13 ] has modified Shu's procedure. He suggested 
that items be ranked in the decreasing order of demand/set-up 
cost ratio and that tho item with smallest ratio be selected 
as skipped item. 
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Andrews and &mmons [ l ] have presented a branch and 
bound algorithm for finding the optimal ordering policy for 
a multi-product system with deterministic demand and no 
backlogging. CO s-fcs were assumed to include a 

separate set-up cost ' for each product ordered but 
witha s fixed savings, -£, each time all products were ordered 
simultaneously. The method, however, does not take care of 
the savings in cost when a few (not all) items are ordered 
at the same point. 

Ignall [ll ] in his paper developed a method to 
determine minimum average cost ordering policies for conti- 
nuously reviewed two product systems with joint set-up costs. 
Markov Renewal Rrog ramming was used to find the region in 
parameter space where a given policy was optimal. 

Herron and Hawley [ 9 ] have considered a multiproduct 
inventory problem where the objective function includes a 
term for storage area cost which is a function of maximum 
expected inventory level G-crald R. Saxton [21 ] has considered 
extensively, product interactions due to storage area. 
Saxton's work consists of developing an analytic method for 
determining an optimum cost policy as a function of order 
quantity for a total cost multi-product inventor model 
where demands are known with certainity. The objective 
function is a sum of storage, replenishment, carrying, 
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obsolescence, area, bin, man-power and equipment costs. A 
computer simulation is included to varify input parameters as 
well as the resultant system cost as provided by the analytic 
model . 

lastly, Arvind R. lele [15] has proposed a model for 
a multiproduct inventory problem with budgetary and store gc 
space limitations, lele structures the situation as a non- 
linear programming problem. The problem is solved by lavidon- 
Fletcher-Powell Variable Metric Method with Golden Search 
technique applied for uni-d irectional minimisation. 

The review of pertinent literature in the field of 
inventory control has uncovered several aspects that merit 
further attention. One of these areas concerns the inter- 
actions which are present in a multi product system. It was, 
therefore, felt that solutions making partial incision into 
this area would be an addition to this field. 

The first part of the present study deals with a 
multi product inventory problem where there is an upper 
limit imposed on total available storage space and items 
are competing for the floor space. The second part of the 
study investigates a problem, where product interactions are 
present due to joint replenishment of two or more items. 



CHAPTER III 


TWO PRO BIB MS : THEIR PORMUIATIOITS AHD SOLUTION 

PEOCEDURES 

In this chapter we consider two multi -item inventory 
problems, and give their solution procedures. The first 
part of the chapter presents a dynamic programming algorithm 
to find out optimal ordering policies for a multi item 
inventory system subject to warehouse capacity constraint. 

In the second, joint replenishment of several items at a 
time is considered and three heuristic procedures are 
developed to solve this prcblem. 

3 . 1 Problem I s 

A statement of the problem is given be lows 

3.1.1 Statement of the problem ? 

A brief description of the decision situation under 
consideration is given b Glows 

Several items are maintained in inventory to meet 
the demands. In satisfying these demands, the systems 
manager finds it necessary to replenish his stock for 
each item periodically. A limited warehouse spaco is 

p 

available for stocking the items. The manager must ddei do 
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when to procure oach item and how much of each item to 
procure so that at no time the inventory level exceeds the 
maximum permissible by the available warehouse space. The 
decision environ mentis composed of item dependent parameters. 
The objective is to dote mine the procurement level and the 
procurement quantity for each item so that the sum of all the 
costs associated with the inventory system under consideration 
is minimized. 

3.1.2 Assumptions ; 

The following simplifying assumptions have been made 
in formulating tho above problems 

1. The demand rate, replenishment rate and lead time for 
each item are known and constant. 

2. Backorders are p emitted. 

3. Item cost, procurement cost, holding cost and back- 
order cost for oach item are known and constant. 

4. Only a limited warehouse space is available. 

5. The maximum inventory level in the warehouse is sum 
of the maximum inventory levels for individual items. 
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.1.3 Notations s 


The fol loving notations arc used to define the problem; 


TSG 


TO, 

TO* 


HC i 

BO, 

PC . 


C i 


r i 

10 , 


b i 


Pi 


m 


i 



= total system cost per period 
= total cost of item * i 1 per period 
= minimum total costof item i per period 
= holding cost of item i per period 
= backorder cost of item i per period 
= procurement cost of item i per period 
= unit cost of item i. 

= demand rate for item i 
= replenishment rate for item i 
= total unit cost of d, items 

= cost of holding one unit of item i in stock 
per period 

= cost of having one unit of item i backlogged 
per period 

= cost of procurement of item i per procurement 
= number of periods per inventory cycle of’ - 
item i 

= lead time for item i 
= Older quantity for item i 
= optimum order quantity for item i 
= procurement level for item i 
= optimum procurement level for item i 
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B i 


¥ 


w i 


= maximum inventory level in an inventory 
cycle of item i 

= maximum backorder level in an inventory cycle 
of item i 

= to tel available warehouse space 
= warehouse space consumed by one unit of 
item i 


3.1.4 Pr ob lem f o rmulation ; 

Since both demand and lead time for all the items are 
assumed deterministic , the resulting inventory process will 
be as indicated in Figure 1. The exhibited geometry of such 
a process will depend on procurement level. Ip, the procure- 
ment quantity, Qp, the demand rate, d ±f the procurement lead 
time, t^p and the replenishment rate, rp. Since these para- 
meters are item dependent, it is evident that the geometry 
will bo item dependent. 
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Math em a ti cal Re lat i onsh i ps 


From Figure 1, it is evident that the number of 
periods per inventory cycle for item i ares 


Also, the following relationships are evidents 


(t u + t 2 i ) (r i - d i ) = + no a -i 


3i 4i' i 


t li + t 2i 


Q i " ^li + t 2i^ ^ r i " d i^ 


t li + t 2i 


t 3± + t 4 . = 


I. + d... t T , - L. 
1 1 lii l 


From equations (2), (3) and (4), we get 


I ± = (1 


) + L i " d i*' fc Li 


The total time per cycle during which inventory is 
hold is (t 2i + t^). This is given by. 


t 2i + t 3i 


h , h 

r. - d ± d . 


From equations (5) and (6), we get 



u- y 


t 2i + i 3i = [Q. (1 - ii ) + I. . + £-) 


We observe the following relationship; 


(7) 


r i 


i + 1 

~ d i d i 


l 


d . 

d. (1 - ~) 
1 r i 


( 8 ) 


From equations (7) and (8), we get, 

*2i + t 3i = [Qq (- 1 - ” ) + L i “ cl i • 


d7 

d- (1 “ 7“) 


(9) 


The total time for which a backorder exists is 
(tii + t^). This is given by i 

B. B. 

*11 + *4i = ?T^d7 + dt" 


( 10 ) 


But, since B i = d.^ . t-r_. - i i , we got: 


*->-• + ~ ^ d i* *li “ L i' 1 ^r, - d, + "d_- ^ ( ' 11 ^ 


li 4i 


i r . 

ill 


From equations (8) and (ll), we get the following 
re lati ons h ip ; 

*li + *4i = ^ d i* *11 ” L i^ ~ d7 

d i ^ - r;> 


( 12 ) 
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Total System Cost ; 


The total cost per period for each item will be 
a summation of the item cost per period, the procurement 
cost per period, the holding cost per period and the back- 
order cost per period, that is; 


TCi = IC i + PCh + HCh + BC ± (13) 

The item cost per period will be the product of the 
item cost per unit and the demand rate in units per period, 


IC, = 0. d. (14) 

ill ' 

The procurement cost per period will be the procure- 
ment cost per procurement divided by the number of periods 
per inventory cycle, 


p. .d. 



(15) 


Holding cost per period will be the product of the 
holding cost per unit per period and the average number of 
units on hand during the period, 


HC i 


1 

T 


b-; ( ^2± + "^-i ) 




m , 


Substituting for 
we get : 


(t 21 + t 3i ), i i 


and in the above 


e quat ion, 



i 


HG ± = 


h i 


c . 

jL < 


d . 


2Q i (1 - jr) 


[«i (i - + h - Vhh 


(16) 


Backorder cost por period will Too tho product of the 
backorder cost per unit per period and the average number of 
units backordered during the period, 


B V = 2 


i_ h (t ii + 


m. 

l 


o r 9 


BGi 


= b i - V 

2Q ± (1 ~ ) 


(17) 


The total cost por period will be a summation of the 
four cost components given by equations (14), (15), (16) 
and (17). Therefore, 


TG i 


< v a ± + 


?l- a l 

Q-i 


h. 


(1 - ^ ) 


d . r\ Id. (L. ■** 

C«i (1 - 5?) + h - Vhd] + — — i 


2 Qi - 57 > 


( 18 ) 


The total system cost will be given by. 
n 

TSG = TG ± 


(19) 
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Warehouse Capacity Constraint ; 

The i-th item in tho inventory system consumes certain 

amount of warehouse space, w^. There exists a certain amount 

of total warehouse capacity, x¥. The maximum accumulation 

of inventory for the i-th item, I_^, will consume 1^ w^ 

units of warehouse space. Therefore, the restriction 
n 

Y~- -^i w i W must not he violated. In the following 
i=l 

section, a dynamic programming algorithm is presented that 
will find the optimal procurement and inventory policies in 
the face of this restriction. 

3.1.5 Solution procedure ; 

The first step in finding the optimal procurement 
and inventory policies for the inventory system, is to 
determine whether tho constraint is active or not. Below, 
we give the procedure for finding tho optimal policies when 
the constraint is active as well as when it is not active. 

1. Optimal Policies when Constraint is Inactive? 

Tho minimum cost procurement level and procurement 
quantity for each item may be found by setting the partial 
derivatives equal to zero and solving tho resulting equa- 
tions. Prom equation (18) we obtain 2 
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°i* d i * 




h i* Q i ^ 


h i 


h i ^ d i #t ii " b i ( d i ,i: Li " L i^ 

+ 57 + d“ 

2Q (1 - -i ) 2Q. (1 - ) 

x r i 1 T i 

( 20 ) 


Taking the partial derivative with respect to Q^, then 
with respect to (d^.t^ - L^) and setting both equal to zero 


give 


^ TC i 


d . 

p r a i , ^ ( 1 ' ? I ) hj(a r t L i - ip 2 

o 2 2 P d. 

Q i 2Q ± " (1 - 

i 


h (a r t L1 - y 

? d. 

2 h a - it > 


(21) 


^(d±.t L . - 1 ± ) 


h l + 


h i tt-ti - L ±> 

d i 

<k (1 - V > 


p (a 1 .t Ii - i t ) _ 

ar~ 

h < x - 57 > 


( 22 ) 


Equation (22) can be written as: 


d • t x • — 1. 
l Lx i 


d i 

h i ^ “ r”^ 
h ± + b. 


( 25 ) 
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Substituting equation (23) into (21), we get 

2p, . d. 




= [ 


x "1 


d. 

(l _ _i ) 
r i 


( i + \ )] 


(24) 


Substituting' equation (24) into equation (23) gives : 


* 

L i 


= W 


- C 


^ • d • ( 1 *p , / -j 

a. j- i ir 

IT ] 

b.. (1 + — ' 


(25) 




J 


Equation (24) and Equation (25) nay now be substituted 
back into equation (20) to give an expression for the minimum 
cost. After several simplying steps, we obtains 


* 

TC. 

l 


= C. 


r d i + [ 


2 (i - % p ± q h ± b ± , 

1 ] (26) 


(h ± + b i ) 


2. Optimal Policies when Constraint is A.ctivc s 


Since I. consumes warehouse snace it is necessary to 
1 

express TC^ for each valueof 1^ over a reasonable range. 

These values make up cost functions for thedynamic programming 
algorithm. 

Equation (5) may be solved for (d^.t^ - Iu) 

giving: 




d. 

Q ± (! - 7“ ) 
1 




(27) 


Substituting equation (5) and equation (27) into equation 
(18) gives: 
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TO. 


p . . d . h . . I . 

,-X d . li 1 . ... 1 1 

°i 9 °i + ^T' £7 


2Q. (1 - 

1 

h [ Q i (1 -^ - y 2 


2Qi (1 - ) 

1 


(28A) 


or 


T0 i 


p . . cl. h. *1 . 

rt ^ , *1 X „ 11 

a i- a i + -q— + ar- 

2Q. (1 ~ •— ) 

i 

d . 

h (l - 57 ) 

+ 5 1 Q ± - 1 ± 


b, .1, 
1 x 


d. 

2Q ± (1 - 


( 28 B ) 


Taking the partial derivative of TQ^ with respect to 
Q. and sotting equal to zero we get; 


^ G i 


Pi* d i 

j- 

Q i 


h ± v 

— T” 37" 

2<V (1 - =r) 

X 


d. 

b. (1 - ™ ) 

X J- i 


b . . I . 
x x 


? d -i 

2Q ± 2 (1 - ~”) 
i 


0 
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This gives; 


Q- 


1 


d ; ? 

2 p. .d.. (1 - + I d (b + h.) 

i 


a . 

(i . -i) 
r ± 


— ] T 


o . 

i 


29) 


From equations (5) and ( 24 ) s we get 


2 Pi-q (i -^) + i ± 2 (h 1 + b ± ) 

L i ~ I i + a i * "^Tji ~ E : ] 


*i* Li 


b. 

i 


(30) 


Equations (29) and (30) give tiio minimum cost and 
tlio minimum cost as a function of and other parameters. 
The minimum cost may be expressed as a function of 1^ and 
other parameters by substituting (29) and (30) into (28A) and 
modifying the last term. This gives; 


TC* 


= C d + + 

i * 1 Q? 

-1 


h. I. £ 

i i 


* d . 

20, (l - zr- ) 


r i 


4“ 


b . ( d . . t_ . 
i v i Li 


r ) 2 


d. 

2Q* (1 - ~ ) 


(31) 


Equations (29), (30) and (31) may now be used to 
develop cost functions needed by the dynamic programming 
algorithm. In the dynamic programming algorithm presented 
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here, each stage corresponds to an item. Let g^ (I w. ) 
represent the total cost when invent oiy level is 1^ and 
lc^ fj_(¥) represent the minimum cost value at the i-th stage 
(i.o. i-th item) when warehouse space available is *¥' 

(¥ is the state variable) , 

The solution proceeds stage-wise with the aid of 
recurrence relations technique, .First, the cost expected 
from the fixot stage (item) if aH warehouse space is 
allocated to it Is determined from f-,(W) = g-^(l-^w^). 

The computations .for f'p(W) are completed and results are 
entered in the first stage of tho solution table. In 
general, for the g'--th stage, following recurrence relation 
is used? 

f-j(W) = Min [g. (I -w ) + f. _ (¥ - I jW . ) ] 

J 0 c I -w . *r V J J J J J J 

(32) 

The results for each stage are entered in the solution 
li'blo. The solution table is then used to find the optimal 
policies for this constrained inventory system. The minimum 
total system cost is searched in the last stage of the solution 
table and corresponding inventory level is noted. Then we 
trace bach through all the stages and inventory level at 
each stage (item) are noted. The optimal values of inventory 
levels so found i.o., I is, along with the equations (29) and 
(30) can be used to find out optimal order policies for this 
inven t o ry sy s tom . 
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5*2 Problem 2 % 

The problem formulation and solution procedures for 
a multi-product inventory system with joint replenishment 
of products are presented in this section. Three heuristic 
procedures have been developed to handle both the cases 
when planning horizon is finite as well as when it is 
infinite . 

5.2.1 Nomenclature ; 

The notations used in this section are given below; 
n = total number of items 

N = total number of set-ups in the planning horizon 

(for heuristic 1) 

= total number of set-ups por year (for heuristics 2 and 5) 
N* = economic number of set-ups per year 

= number of set-ups for i-th item in the planning 
horizon (for heuristic l) 

_ number of sot -ups per year for i-th item (for heuristics 
2 and 5) 

kt* = economic number of set-ups for item i in the planning 

i 

horizon (for heuristic l) 

= economic number of set-ups per year for item i 
(for heuristics 2 and 5) 

k ± = relative ordering frequency- for item i (k ± = 
k* = economic value of k. 

T 1 
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relativs ordering frequency for item j when item i is 
ordered at each set-up 
corrected value of k.. 

relatire ordering frequency of item i at the r-th 
iteration 

corrected value of k^ r 
planning horizon 

number of periods in the planning horizon 
number of time units per year 

time duration in an inventory cycle for which a. unit 

of item i is backlogged 

order quantity for item 1 i’ 

economic order quantity for item i 

annual demand for item i 

demand rate for item i 

fixed cost of talcing a sot-up 

variable cost of taking a set-up for item i 

cost of holding one unit of item i for unit time 

cost of having backlogged one unit of item i for 

unit t ime 

warehouse space ^occupied by one unit of item i 
total available warehouse space 
cost of occupying one unit of warehouse space 
during the planning horizon (for heuristic l) 
annual cost of occupying one unit of warehouse 
space (for heuristics 2 and 3) 



YG i 


TC 



variable cost for item i in the planning horizon 
(for heuristic 1) 

annual variable cost for item i (for heuristics 2 
and 3 ) 

total cost of the inventoiy system in the planning 
horizon (for heuristic l) 

annual total cost of the inventory system (for 
heuristics 2 and 3) 

annua. 1 total cost of the inventory system when item 
i is ordered at each set-up 


3.2.2 Case of Finite Planning Horizon s 
1. Problem Statement ; 

There aro n items stocked in a warehouse. The various 
costs, namely, 

(i) Fixed and variable set-up costs, 

(ii) Holding cost, and 

(iii) Backorder cost 

arc known for each item. Decisions as to how many orders arc 
to bo placed, when to place them and how much to order each 
time aro to bo made. Those decisions should minimize the 
total cost incurred over the finite planning horizon, T. 


2. Assumptions ; 


The model developed, is based on a sot of assumptions. 
The purpose of these assumptions is to simplify tho analysis 
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of the problem. These assumptions are described below; 

1. There is a family of n products which can be ordered 
together. 

2. -An item is replenished at equal time intervals. 

3. Backorders are permitted. 

4. Replomis hment is instantaneous for each item. 

5. There is only a limited storage space available for 
stock holding. However, additional space can be 
rented at a constant rate. If there is some extra 
space it can be rented out at the same rate. 

6. The planning horizon is finite and known. 

7. Demand rate for each item is given and is constant 
ove r t ime . 

8. The demand for each item within the horizon period is a 
quantity greater than zero. 

9. There will be at least one set-up for each item within 
the horiz cm period . 

10. The planning horiz on consists of a finite number of 
periods and a set-up can be taken up only at the 
commencement of a period . 

11. An order will be placed for each item at the beginning 
of the first period. 
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13. The numb cr of inventory cycles for each item, within 

the planning horizon, is an integer. 

3. P rob lem Po rm ulat ion i 

The inventory process resulting from the problem as 
stated above is shown in Figure 2. 



Fig. 2: Inventoxy Process for Problem 2 (Finite Horizon) 

Since the i-th item is replenished in equally 
spaced set-ups in the horizon period, the variable cost for 
this item consists ofs 

(i) variable set-up cost of N ± set-ups? 

(ii) total stock holding cost in the plaming horizon? 

(iii) total backorder cost in the planning horizon? and 

(iv) cost due to ware house space occupied. 


Thus, we can write, 
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VG i = l^.s. + i (Q i -t ± .d i )(|- - t ± ) ^ ± .h. ± 


+ £ t?od i .IL.b i + (Q i ~t i .d i ) Wj^.R 


(33) 


Tho total cost of the inventory system can ho written as, 


n n 

10 = K3 + N.s. + £ y~ (Q.-t..d.) x 

iS 1 x i5T 1 1 1 


n 


- V v*i + * 


n 

+ R ( T_ ( CL -t . . d . ) w. - ¥ ) 
i=l 111 i 


(34) 


where the order quantity, Q^, is given by 

Td. 

«i = Iff (35) 

To find the optimum value of t- which corresponds 
to iainimum variable cost, substitute the value of in 
equation (33) and equate the partial derivative of VCh_ 
with respect to t^ to zero . This gives , 

\r..R + T.h. 

t. = -i (36) 

N i (h i +b i ) 

Substituting this value of t ^ in equations (33) and (34), 
and replacing Q i by Td^/ih , we get, 
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VC i 


F i- S i + ’2N. 


T 2 d. h.,b. 

1 2 . X 


h.+b. 

i i 


+ R 


2N. (h. + b. ) 

i l i' 


and 


n 


TC = 


Id. h..b. 

us + T" 


i=l 


i i 


w .d (2T.b. -w..H) 

+ R “X Th~ ~b ~) ~~ 3 - H.V 


(37) 


(38) 


4 . Solution Pro oo dure ; 

Heuristic Algorithm Is 

The solution procedure consists of two- parts. 

The first part consists of finding the values of for each 
item. In the second part, these values of 1 are plotted 
and all the points at which one or more items are ordered 
arc counted. This gives the total number of set-ups, in 
tho planning horizon. 

The economic ordering frequency, for the i-th 

item corresponds to tho minimum variable cost for the i-th 
item. The vaiuc-s of it for i-th item can, therefore-, be 
obtained by taking the first derivative of equation (33) with 
respect to E ± and equating it to zero. This gives, 
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IT, 


T .dj_.hj.-bj + R.w i .cI i (2T.b i - Wp.R) ^ (35) 

2s i (ll i + V 


This value of Kp is then corrected to ensure an 

integer number of invent oiy cycles in tho planning horizon. 

A stop by stop algorithm of the proposed heuristic 

is given belowj 

Par t 1 % This part consists of obtaining the values of Up. 

Step 1; Prom the given value of P find out the possible 
values which N- can take. Since the number of 
set-ups, lh, for the i-th item should be on integer, 
tho possible values which Fp can take will be tho 
factors of P. let these bo denoted by a 2* * ° a p ! * * 
Thus, for each item i the number of set-ups for that 
item will be given by one of the a p ’s. To illustrate 
this, let P = 12, then the values of a p will be 
a^ = 1 , &2 ~ 2 5 a 3 = ^ > a 4 = 4 , = 6 and a^ = 12 . 

So Fp can take any of the values 1,2, 3, 4, 6 and 12. 

Step 2: Find out the va.lue of Fp for i-th item from 
equation (39)* 

Step 3: If this value of Up is less than one then put Fp = 1 
and go to Step 9, otherwise continue. 

Stop 4s If the value of F ± calculated in Step 2 is greater 
than P , then put F* = P (This is done because there 
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cannot be more then P set-ups in the planning 
horizon) and go to Step 9 otherwi se continue. 

Step 5; If the -value of N ± is equal to one of the factors 
of P, say a^, then put = a^ and go to Step 9 5 
otherwise continue. 


Stop 6 s 


If the value of N. lies between a. and a-,-, thon 

1 J J'T* — 

t 

put 1'R = a. and compute the value of variable 
cost VC il s for i-th item by the following 
formulas 


T.d. h . . b . 
VC il = F r s i + 2ir- h~+b~ 


w..d . (2T.b . - w .R) 

- L -_ rb. - _.4- - jL - - 

(h i +b i ) 


Step 7: 


Take the next higher value of a i.e. a^ + ^ and lot 
NV = a -j + y Again compute the variable cost 
with this value of L . let it be VC^ 2 » 




+ 


T ? .d ± hj .b^ 
211 " h.+b. 

l i l 


+ R 


w i .d jL (2T.b i - w ± .R) 

2n" VV 


If the value- of VC ±1 is less than or equal to 
VC ±2 then put IP = 1 $' ± and go to Stop 9, 


Step 8: 



33 


otherwise put 1^ = F*\ 

Step 9; Repeat Steps 2 through 8 for all items. If the 

values oi 17^ for all the items have "been calculated 
then go to Step 10. 

1 P This part of the algorithm finds out the total 

number of set-ups , N, th ©total cost of the s; r stem 
and the economic order quantity for each item. 


Stop 10: Using the above calculated values of find 

out R. for the i-th item whore R. = p/R.. 
i i ' i 


Stop 11: An order for evoiy item is placed at the first 
period and successive orders are placed after - 
every R^ periods. Rind out tho periods at wliich 
an order is placed for tho i-th item, let those bo 


denoted by p^, p i2y . p^,... p.^ (where the 

maximum value of p^ can be P and minimum value 
can be l). 


Stop 12: Repeat Stop 10 and Step 11 for all tho items and 
find out the values of p . 1 s for all the items. 

Stop 13: Count those periods at which one or more itoms 

are ordered. Tho total number of such periods give 
the tote.1 number of set-up" y IT. 

Stop 14: Obtain tho economic order quantity for i-th item 
fran the following fo rmula: 



Stop 15 J Obtain the total costof the system by substituting 
values of IT snd Iff , in the equation (38) . 

3.2.3 Case of Infinite Planning Horizon ; 

In this section two procedures are described 
which deal xii tb a joint replenishment problem of a multi - 
item invent ory system when the planning horizon is infinite. 
The statement of the problem is given below? 

1. Problem Statement ; 

In the previous model we treated the planning 
horizon as finite. The same problem is dealt here with an 
infinite planning horizon. The decisions are to be made 
such that the total cost per year is minimum. 

2. Assumptions 1 

The first five assumptions made in the previous 
model hold true- for this case also. Following additional 
assumptions have been made? 

1. The set-ups are taken at equal time intervals. 

2. Planning horizon is either infinite or it is 
so long compa.ro d to the cjrcle time that it 
can betaken as infinite. 



3. Atleast one item is ordered at every sot-up. 

4* Total number of sot-ups per year need not be an 
integer . 

5. All items arc ordered at the first set-up. 

6. Annual demand for each item is a quantity greater 
than zero . 

3. Problem Formulation r 

The inventory process resulting from the problem 
statement given above is shown in Figure 3- 



Fig. 3 : Inventory Process for Problem 2. (Infinite 
horizon). 

If the i-th item is replenished at every k^-th. 
set-ups then there will be H/k^ set-ups per year for this 
item. Hero, k ± must be a positive integer and is called as 
relative ordering frequency for the i-th item. The annual 
variable coot VC^ for the i-th item is given by: 
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VG ± 


w t . ,D. 

— - s . + i ( C _ x l 

k ± i 2 (f *i 5 


mk. 

> <t - *-■ > 


IT 


x' k ± 




2 ^ t . . D . 

+ * h Vk~ h + E ' w i 


m 


(40) 


Ths total annual cost forth© inventory system is given by 

a ^ n 


TC 


NS + 


1=1 i 


w j--l t . . D . mk . 

s_. + 4- (Q, - 


k_. i 


i=l 


m 


li 


TSS 


n 


2 D i- H 


1 S i ’ 2 111 V ^7 “1 + R[ w i (Q i~ _ lr ;) " W l 


where the order quantity, Q. for item i is given by. 


Qi = 


D.k . 

x x 

IT 


(41) 


(42) 


In the expression for VC^, substitute the value of 
from equation (42) and equate the partial derivative of 
VCh with respect to t^ to zero. This gi'ves optimum value 
of t^ as: 


fik.v . + n.k. .h. 
t„. = —i— i — - 

1 N(h ± + b ± ) 


(43) 


Substituting this value of t. in equations (40) and (41) 


and replacing Q. by D-.k^/lT. , we get: 


l 

ITs 


VC ± 


D. ,m 
_i + 1— 

k ± + 2 IT 


h. .b. 

_i 1 v 

h ± +b i i 


w. .D. .k. ( 2n.b . - R.w. ) 

+ R kt 1ST A ± ) 


( 44 ) 
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and , 


TG 


, lf , JL r IIs i D-.u h..b. 

1,10 + 2 L ir + *~orT" • -u " ~~ v ^ k + 

i=l 4F h i +b i i 


+ „ w l- B i- k i ( 2ia - b i - Bw ± ) 

ii: I (h ± + b.l J - 


R.W 


(45) 


G o lu't ic n I to c p dura ; 


Heuristic Algorithm 2; 


In the above p rob Ion , two policy decisions arc 
to be t ;ikou : ths, total number of sot~upts per year N and 
the number of sot -ups for individual items N/ k ± . For a given 
value oi J, the tern NS in equation (45) is constant. The 
valu s, of k| for the i-th item can then bo found by differen- 
tiating equ,' lion (44) partially with respect to k^ and putting 
it equal to zero. This gives s 


N [ 


28 r m (h ± + b ± ) 


I>. #n \h. .b. 4- R# ¥ . . D. (2mb- - fLw. ) 

1 i- X 11 1 1 


P (46) 


or. 


if 

k~ 

i 


Ih . ,h. .b. + R.w. .D. (2mb.-R.w. ) ± 


l l i i v i i' 

2a i .n (h ± + b ± ) 


] 2 (47) 


We have assumed that at least one item 'will be ordered 
at each set-up. This means at least one of the k^'s will be 
equal to 1, The procedure adopted here can now be explained 
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as follows; First find out the values of N/k^ for each item 
from equation (47). Since one of the k^s should he 1, 
therefore, start with, the first item and put the correspond- 
ing value of equal to 1. This gives a particular value 
of E. From this value of E find out tho values of k ± for all 
other items. Calculate the total cost of the system with 
those values of IT and k^'s. Next for the second item put 
k^ = 1 and repeat the whole process. The above procedure is 
followed for all items and the one which gives minimum total 
cost is selected. 

A step by stop algorithm for the abovs procedure is 
given below; 

Step 1: Use the expression (47) to find out the value of 

F/k^ for tho i-th item, let this value of N/k^ be 
called as N - . Repeat this process to determine 
values of for all items. _ 

Stop 2; Atloast one of the items is to bo ordered at all the 
setups or in other words it can be said that atleast 
one of the k^*s should be one. let it be for the 
i-th item i.o. , k^ = 1. So the value of N (i.e. total 
number of sot-ups) will bo tho same as N^. 

Stop 3; If for any other item j ( j 4 i) tho value of IT . is 

J 

greater than or equal to the value of IT calculated in 
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Step 4 


Step 5 : 


Step 6 


Step 2 then for that item put k c . =1 and go to 

J a 

Step 9, otherwise continue. 

Calculate the value of k.. for those items only for 
which Nj is less than Ep The value of k.^ for such 
an item can he calculated from the relationship 
given below; 

N. 

k . . = —■ 

ji N . 

J 

If for any item j, the value of k.. is an integer 

J 1 

thon put k° = k n ... and go to Step 9, otherwise 

J -h J -h 

continue. 


If kp is a fraction then an integer value of it is 
to bo taken, which gives the minimum variable cost for 
the j-th item. In such a case k.. can take either 

j “k 

lower integer value, soy k .-p or the next higher 
integer value, say k^ i2 . ^isst take the lower integer 
value and calculate the variable cost, soy VC 

t. % ' L * 

for the j-th item with the value of E calculated in 
Stop 2. The value of VC ^ can bo calculated by 
using the following equation; 


E. . s . D . .m 


TC ji 


— -Jtl + tl 




h ..b . 

1 k 

2^ * h +b ‘ k jil 


R.v ..i)..k. i:L (2m.bj - R.w,.) 


2m.h i (h. + b.) 

J J 
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Step 7: 


k^2 ? and again 


Take the next higher integer value, 
calculate the variable cost, say VC.g, for the 
same item, using the following equations 


VC 


32 


N i <s .i 

k d±2 


D . .m 
-4.L 




k . . 


2T ± * h j+b j ‘ j'i2 


R.w ..D . .k . . ? (2m b. ~ R. w • ) 

+ « J rl a±rA A u — . 

2m. I. (h. + b.) 

1 J J 

Step 8s Compare the values of VC.-^ and VC If VC.^ is 

J c 2 

greater than or equal to VC.^ "then pat 

c f 

and go to Step 9, otherwise, x put k^ = k^. 

Stop 9: Repeat Steps 3 through 8 for all the va3r.es of 3 
(i.o., for the value of IJ calculated in Step 2 
find out the values of k^. for a3JL tho items). 

c 

Stop 10: Calculate tho total cost for the set of k^'s and 

R values calculated above. The total cost, TC^, can 
be computed using the following equation; 


TC. 


n 

N.S + S' 

1 J=1 


C 


N.s, D..m h . . b . 

^ .1 + _3 ..3.-JL . k c . . 

c + 2N. * hi+bT 31 


k.. 

31 


3 


R.w ..D-.k'J. ( 2mb . - R.w.) 

3 — 3 — li 3 R>w 

2m N ± (h • + b j ) 


Step 11; Ropcat tho Steps 3 through 10 for i - 1,2,..., n 
(i.o. putting k ± = 1 for the next item calculate 
the values of N and k^’s and corrosponding total cost). 
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Step 12: Select that value of II and set of k c .'s for which 
the total cost TCh is minimum. let the selected 
values of k^ . *s he represented as k. 1 s for all 
items . 


Step 13 % 


Calculate the economic order quantity for all the 
items using the following expression: 


Heuristic Procedure 3: 

The method described in this section is an iterative 
procedure. To start with, a set of values of are arbitrarily 
selected and the value of H is determined. The corresponding 
value of the objective function is then calculated. The value 
of the objective function is improved at each iteration. If 
for any two consecutive iterations, value of the objective 
function remains tho samo, then the procedure is stopped. 

The logic behind the proposed procedure is- as follows: 

For each sot of values of k^'s there exists a local 
minimum. The value of H which corresponds to this minimum 
can bo obtained by taking partial derivative of TC, in 
expression (45) with respect to N and equating it to zero. 

This gives: 
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N 


2 , 

m > 


n 


3, . 


h. .b. 
x 1 


± 

T 


i=l h i +l3 i 


n D . . w . k . ( 2mb . -R .w ) 
k +Rj-i 11 1 1 

i=I 


( h • + b . ) 
1 x 


l— 


2m 


( S Hr 


n 


i=l 


fi. 

V 


J 


(48) 

Again, a, local minimum exists for every value of IT and corres- 
ponding values of k ± 's for all items, for a given value of IT, 
the term ITS in the expression for total cost is a constant. 

The optimal value of k ± for item i can, therefore, bo 
determined by tailing the partial derivative of variable cost, 
VCh, with respect to k^ and equating it to zero. This gives s 


k. 


IT [ 


2sj_.m (bn + b^) 


D. .Tsr h. .b. + R.W..D. (2mb. - R.w.) 

X 11 111 1 


i 

T2" 


The economic number of set-uxps., ST, and re la- l 

tivc ordering frequency k^ for the i-th item can be deter- 
mined with the help of expressions (48) and (49). The value 
of IT depends on a prior knowledge of k^ for each item, which 
in turn depends on IT. An iterative procedure is, therefore, 
adopted for determining IT and k^ for individual items. 

For the r-th iteration in the iterative procedure, 

we . 


(a) use the k. values obtained in (r-l)-th 

-L y Jl “"-L, 

iteration when IT was equal to IT and 
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(b) Obtain K r and r £or each item. If the value 
of k i,r is not 311 integsr, then take the lower 
integer as value of k. and calculate the 

1 $> Jl 

value of variable cost for i-th item, let it be 

VCh , 1 * take the higher integer as value 

ot r and calculate the variable cost, let it 
# ^ 

bG VO i,2* If VC i,l is lGSS than Y0 ± o then the 

lower integer is taken as -the value of k. 

i,r J 

otherwise, the higher integer is taken as k. 

3- f T 

The iterative procedure will be terminated in r-th 
iteration if k ±?r = ^ x for all items. 

The following steps arc involved in the proposed 
iterative procedure; 

Stop 1; For the first iteration put r = 1 and select arbitrary 
values of k ir for all items. In practice, it is found 
that values of 3sa r for i = 1,2, ...n, provide a faster 
convergence. (If the arbitrarily selected values of 
k. ,’s happen to be near the optimal values of k. * s 

J-.J1 XT 

then the convergence with those values of k^ r f s will 
be faster. ) 

Step 2; Calculate the total number of set-ups IT from the 

formula given below, lot it be represented as ET . 

Thus, we gets 
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Step 3; 
Step 4 : 


Step 5 


Stop 6 

Step 7 


o n h . . b . 

~~ Tj -i 1- ir . t) <c- i 7 Tjr 

i=i W ir ‘ i=i rE~^TT7) 


n D n . (2mb^ - R. w^) i 


1 

1 "2T 


n 


s . 

i 


2m (S + yj ) 

i =1 ‘"ir 


(50) 


Set r = r+1 and i = 0 , 


Set i = i J -l. Use the value of N r calculated in Step 2 
(i.e. in the previous iteration) to find out the new 
value of k^ r « The value of k^ r can be calculated from 
the formula given below % 


k 


2 b ± ) 


k. 


N 


r-1 L 1 _2 


D. • «b .4- R.w - . D. ( 2mh_. -R.~w . ) 

1 JL JL 1 jl * j.. i. 


] T (51) 


If the value of k^ r calculated in Step 4 is less than 
or equal to 1, then put k° r = 1 (This is done because 
no sot-up can be taker at intermediate points) and 
go to Stop 10. If the value of y is greater than 1 
then go to Stop 6, 

; If k? r is an integer number then put k?^ = k^ and go 
to StcpilO, otherwise ? continue. 


Take the lower integer as value of k ir , 


let it b e 


repro sente d as k^ r . Calculate the variable cost for 
i-th item from the formula gi-ven below, let it be 
VO^. Then W ±1 is gi-ven by: 
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vc., 

1.1 


h-l 


^ir 


D,. . m h . . b . 

_i 1 k * 

2 V-l' h i +b i ‘ lr 


R * w i- D i* k ir “ R.w ± ) 


2m. H 


r-1 


(h, + b, ) 


Step 8; Take the next higher integer as value of ^j_ r * Bet 


t? 


it be represented as k. . 

lr 

cost with this value of k. 
V 0^2 is given by: 


Calculate the varis.ble 
. let it be VC^* Then 


N 


VC 


r~l ' 


i2 


D. .m 

+ ^ 


ir 


r-1 


h..b. 
i i 

h.+b. 
i i 


k!’ 

ir 


R.w i .D i .k il , (2mb i - R.w i ) 

2m. iV^ (h ± + b ± ) 

Step 9: Compare the value of VC i1 with VC.^ • If T ^il ^ 

C f 

less than VC i2 » then put k^ r = k^ r and go to 
Step 10 s otherwise , put k? = k" r « 

Step 10: If the value of i isless than total number of 

items n, then go to Step 4, otherwise, continue. 


Step 11: If k? = k. , f 0 r i =1,2,..., n (i.e. if the two 

1 JL 1 J 1 *"* — Li 

consecutive values of k ir are same for all the items) 
then go to Step 12, otherwise , put k ir = k?^ for 


1 — 1.2, ... 


n and go to Step 2. 



Step 12: The economic number of set-ups per year tf is equal 

"k° ^r-1* out the economic order quantity for 

each item, using the following formulas 

D . . k . 


Step 13 


Qt = 


H 


* 


Determine the number of set-ups for each item from 
the following formula,? 


H* = 




ir 


i = 1,2, 


Step 14s Substitute the values of II and L fo: 

into equation (45) and calculate the total cost of 
the system. 



CHAPTER IV 


MB THODO IOCS' TESTING AND SCOPE FOR FUTURE WORK 

In this chapter results pertaining to the performance 
of the proposed models on a few test problems are presented. 
At the end of the chapter, possible avenues for further 
work are sugg est e d . 

4.1 Performance of Pvnamic Programming Procedure : 

Nine problems have been solved using discrete dynamic 
programming computer package and for each problem optimum 
ordering policies have been determined. The problem size 

is varied 'from two items to four items. For each problem 

size a set of three problans have been solved. The average 
computation time for these problems on the IBM 7044 Computer 
system are given in Table 4.1. For tho sake of illustra- 
tion, the complete results of only one problem, for which 
the data is given in Table 4 . 2 , are presented. For this 
problem, the total available warehouse space is taken as 
16 units. The results for the optimum ordering policies 
(i.e. , procurement quantity and procurement level for each 
item) are given in Table 4.3. 

The same problem has also been solved using continuous 
dynamic programming procedure. The computer package by 
Mize [14] was utilized for this purpose. The optimum 
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ordering policies for the problem, as obtained by the conti- 
nuous dynamic programming procedure are given in Table 4.4. 

It is found that the discrete dynamic programming 
procedure takes a computational time of 2.2 seconds and. 
gives an objective function valueof 284.689 units. The 
continuous dynamic programming results in a better objective 
function value of 278.736 units but it takes 17.605 seconds 
of computational time. Therefore, it is inferred;-- that 
continuous dynamic programming approach yields better 
results, however, at a higher computational time. 

4.2 Performances of Heuristic Algorithms. ; 

The performance of the three heuristics discussed in 
Chapter 3 have been compared by solving about 150 problems 
of various sizes. The problem size ism varied from 6 items 
to 40 items. For each problem size a set of problans have 
been generated randomly. Computer packages have been developed 
in FORTRAN IV for the IBM 7044/1401 Computer system, for 
testing the proposed heuristics. 

The performances of the heuristics are evaluated based 
on cost index and average computational time for a given 
problem sizo. The details of the two measures of evaluation 
arc given in the following paragraphs* 
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Cost Index; 


For a given problem size, n , let TC and 10. be 

en m 

the total system cost using economic order quantity model 
and heuristic approach ' i* respectively. Then the cost 
index, CI^, is determined using the following relationship: 


Cl 


in 


TO. 

in 

tc 

en 


wher-c, i - index of heuristic (i = 1,2,3) 

n = index of problem size ( n = 6,7,..., 40) 


For a given problem the heuristic which gives least value 
of the cost index is -taken as the ' bostt heuristic. 


The threo heuristics have boon compared for a set of 
45 problems chosen rand only. Heuristic one gives 'tho total 
system cost for a given planning horizon, whereas, heuristics 
two and three give annual total system cost. To compare 
heuristic one with heuristics two and throe, the planning 
horizon for heuristic one is taken as one year. The results 
are presented in Table 4.5. The results indicate that 
heuristic three yields least cost index values in most of 
the cases (42 out of 45). In three cases the cost- index 
values obtained by heuristic two are least and in 11 cases 
these are very near to those obtained by heuristic three. 

The higher cost index values for heuristic one can be 
attributed to the following reasons which accrue frem the 
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tvro constraints involved in this heuristic; 

there is an upper limit on the total number of set-ups 
per year, 

2. The total number of inventory cycles for each item 

has to be an integer. 

Further, examination of Tabic 4.5 indicates that the 
cost index values obtained for each of the heuristics is 
loss than 1. This implies that it i s economical to replenish 
several items jointly rather than ordering them separately. 

Average Computational Time ; 

For all tho generated problems, the average computa- 
tional time is determined for each problem size. Average 
computational time values for different problem sizes are 
presented in Table 4.6. It can be seen from the table that 
for small si 2 e pr<±> lems (upto 7 items) heuristic two takes 
tho least computational time. Heuristic three closely follows 
heuristic two. For large size problems (from 8 items to 
40 items) heuristic three takes least time. A possible 
explanation for this behaviour of heuristics is given below; 

Heuristic two generates multiple solutions (equal to 
the number of items) for a given problem. Out of these 
solutions, the minimum objective function value solution is 
selected. Thus the computational time of heuristic two 
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incroasos as the problem size increases. Heuristic three, 
on the other .hand, is an iterative procedure, and as will 
"be shown latter, the number of iterations do not depend on 
the problem size (the number of iterations, in general, are 
between 6 to 10). Therefore, for largo size problems 
heuristic throe takes much less time compared to heuristic 
two . 

Based on the two measures of evaluation, it can be 
concluded that heuristic thiee is the most efficient 
heuristic . 


Heuristic three, being an iterative procedure, in- 
further tested to evaluate for its convergence. This is 
done on the basis of number of iterations it takes to solve 
a problem of given size and the percentage improvement in 
the total system cost. 


Bora problem of size 'n', let TC ln and TC Ln represent 
the total system cost at the first iteration and at the last 
iteration respectively. Then the p ercentagenximprovement, 
PI n , is given by the following relationships 


PI 


n 


TC i* 


.TC 


rn 


x 100 


Table 4.7 gives the values of PI n for ^52 Problems 
ranging in size from 6 items to 4^ items. The second 
column of the table gives the number of prol|l^m|.. ,§(317^^ 

ce ii-; 

Acs, 




52 


for the problem size under consideration. The third and 
fourth column of the table give the minimum andmaximum 
values ...iOf the number of iterations and percentage improve- 
ments respectively. As can be seen fran the table, the 
number of iterations in most of the cases are less than 
10 (only In 4 cases the number of iterations exceeded 10). 

Thus, it can be concluded that number o f.i ;ite rations do not 
depend on the problem size. Moreover, considerable improve- 
ment in the ob jective function value can be made in a few 
iterations. To illustrate it, for an 11 item problem 
percentage improvement of 88^13 is obtained in 9 iterations. 

In the previous paragraphs -we have compared the three 
heuristics when the planning hbrizon is infinite. Now we 
shall evaluate the performance of heuristic one which considers 
finite planning horizon. The evaluation of this heuristic is 
done using the following procedure. 

The total system cost is calculated using the heuris- 
tic. The total system cost is also calculated when all the 
items are ordered at each period. The total system cost 
values using both the procedures are compared. 

For a problem of size, n, let TC^ represent the 
total system cost when all the items are ordered at each 
period and let TC ln represent the total system cost calcu- 
lated by heuristic one. The percentage cost saving, FS n , 
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obtained by using.! heuristic one is defined as; 

TC - TC, 

PS n = — — ~ — x 100 

An 

where n = index for problem size. 

Table 4.8 gives the values of PS n for 13 1 problems 
for the problem sizes ranging from 6 items to 40 items. 

Prom the table, i t is clear that the use of heuristic one 
as compared to the procedure of ordering each item at 
every period in the p laming horizon results in a saving. 

For the selected set of problems, savings upto 97. 08- -percent 
were obtained by the use of heuristic one. 

4.3 Scope for Future Work ; 

In developing the methodologies presented in this 
dissertation, a few avenues for further work emerged. These 
are given below: 

1. In the dynamic programming procedure used it is 
assumed that the demand and lead time are d ete minis tic . It 
will be worthwhile to extend it to the case where demand 
and/or lead time are stochastic. 

2. In heuristics two and three it is assumed that there 
is no limit on the total number of set-ups per year. It 
would be interesting to analyse the problem when there is 
an upper limit on the total number of set-ups in a year. 
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3. In the three heuristic procedures developed in this 

work, it is assumed that product interactions only due to 
joint replenishment of the products are present. It would 
he worthwhile to consider the cases where additional product 
interactions due to warehouse capacity restriction and/or 
budgetary restriction are also present. 



Table 4.1s 


55 


Average Computational Time for Problems 
Solved by Discrete Dynamic Programming 
Procedure 


S.No. 

Problem Size 

Average Computa- 
tional Time (Sec.) 

1 

2 

2.333 

2 

3 

3.075 

3 

4 

3.992 


Table 4.2: Data Used for a Three Item Inventory Problem 
Solved by Both Discrete and Continuous 
Dynamic Programming Procedures. 


Item 

No. 

Demand 

Rate 

Item 

Cost 

Holding 

Cost 

Backorder 

Cost 

Pro cu- 
rem ent 
Cost 

lead 

Time 

Space consum- 
ed by one unit 

1 

6 

30.88 

0.30 

0.30 

18.30 

2 

4 

2 

4 

18.33 

0.24 

0.17 

17.50 

4 

2 

3 

1 

12.00 

0.12 

0.25 

15.50 

1 

1 
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Table 4.3; Optimal Ordering Policies Obtained by 
Discrete Dynamic Programming Procedure 
for the Problem Griven in Table 4.2 


Item 

No. 

Procure- 
ment level 

Procurement 

Quantity 

Warehouse space 
allotted 

1 

-14.08 

27.10 

4 

2 

-10.06 

29.07 

6 

3 

- 6.31 

13.31 

6 


Total Cost; 284.689 units 


Computer Time Taken 5 2.2 seconds 


Table 4.4: Optimal u r dering Policies Obtained by 

Continuous %namic Programming Procedure 
for the Problem given in Table 4. 2. 


Item 

No. 

Procurement 

level 

Procurement 

Quantity 

Warehouse space 
allotted 

1 

-12.91 

27.26 

9.400 

2 

-11.20 

28.80 

3.200 

3 

- 8.80 

11.28 

*•932 


Total Cost: 278.736 units 

Computer Time Taken: 17.605 seconds 


Table 4.5: Performance Comparison of Three Heuristics 
Based on Cost Index Values. 


SI. 

No. 

Pr oblem 
Size 


Cost 

Index Values 


Heuristic 

1 

Heuristic 2 

•heuristic 3 

1 

6 

0.482 


0.447 

0,331 

2 

7 

0.637 


0.610 

0.509 

3 

7 

0.904 


0.889 

0.887 

4 

8 

0.704 


0.636 

0.634 

5 

9 

0.723 


0.637 

0 . 636 

6 

10 

0.309 


0.288 

0.202 

7 

11 

0.482 


0.477 

0.420 

8 

12 

0.485 


0.478 

0.421 

9 

13 

0.421 


0.320 

0.249 

10 

14 

0.246 


0.230 

0. 216 

11 

15 

0.320 


0.294 

0.217 

12 

16 

0.591 


0.594 

0.592 

13 

16 

0.960 


0.V68 

0.944 ‘ 

14 

17 

0.211 


0.187 

0.163 

15 

17 

0.318 


0.324 

0.283 

16 

18 

0.278 


0.249 

0.151 

17 

19 

0.591 


0.580 

0,555 

18 

19 

0.982 


0.953 

0.922 

19 

20 

0.527 


0.527 

0.486 

20 

21 

0.374 


0.363 

0.350 

21 

22 

0.723 


0.538 

0.487 

22 

23 

0.641 


0.6 23 

0.6 23 

23 

24 

0.623 


0.549 

0.556 

24 

25 

0.583 


0.548 

0.554 

25 

26 

0.712 


0.698 

0.488 

26 

27 

0.281 


0.299 

0.302 

27 

28 

0.136 


0.117 

0.115 

28 

29 

0.163 


0.146 

0.137 

29 

30 

0.124 


0.114 

0.112 

30 

31 

0.143 


0.149 

0.129 

31 

32 

0.208 


0.187 

0.103 

32 

33 

0.189 


0.153 

0.153 

33 

34 

0.168 


0.131 

0.118 

34 

35 

0.178 


0.246 

0.135 

35 

35 

0.507 


0.532 

0.482 

36 

36 

0.126 


0.116 

0.112 

37 

36 

0.667 


0.675 

0.632 

38 

37 

0.292 


0.284 

0 . 169 

39 

37 

0.862 


0.831 

0.782 

40 

38 

0.224 


0,164 

0.103 

41 

38 

0.703 


0.675 

0.621 

42 

39 

0.152 


0.109 

0.099 

43 

44 

39 

40 

0.578 

0.173 


0.527 

0.162 

0.483 

0.160 

45 

40- 

0.806 


0.713 

0.629 
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Table 4.6s Performance Evaluation of the Three Heuristics 
Based on Average Computational Time. 


SI. 

No. 

Problem 

Size 

Average Computational Time (Sec.J 

Heuristic 1 

Heuristic 2 

Heuristic 3 

1 

6 

0.378 

0.267 

0.270 

2 

7 

0.420 

0.291 

0.295 

3 

8 

0.480 

0.368 

0.342 

4 

9 

0.503 

0.434 

0.320 

5 

10 

0.984 

0.483 

0.429 

6 

11 

0.729 

0.617 

0.392 

7 

12 

0.820 

0.723 

0.515 

8 

13 

0.793 

0.936 

0.581 

9 

14 

0.840 

0.927 

0.601 

10 

15 

1.082 

1.038 

0.642 

11 

16 

0.960 

1.026 

0.687 

12 

17 

1.020 

1.015 

0.730 

13 

18 

1.082 

1.061 

0.720 

14 

19 

1.146 

1.127 

0.761 

15 

20 

1.208 

1.201 

0.804 

16 

21 

1.403 

1.333 

0.785 

17 

22 

1.308 

1.293 

0.770 

18 

23 

1.382 

1.321 

0.920 

19 

24 

1.458 

1.431 

0.963 

20 

25 

1.489 

1.471 

0.916 

21 

26 

1.621 

1.599 

0.923 

22 

27 

1.672 

1.660 

0.978 

23 

28 

1.998 

2.016 

1.023 

24 

29 

2.123 

2.088 

1.042 

25 

30 

2.004 

1.950 

0.857 

26 

31 

2.087 

2.021 

0.885 

27 

32 

2.135 

2.080 

0.913 

28 

33 

2.263 

2.145 

0.947 

29 

34 

2.371 

2. 212 

0.963 

30 

35 

3.398 

3.375 

1 . 021 

y w 

31 

36 

2.367 

2.341 

1 .028 


37 

2.428 

2.417 

1.056 

J c- 

33 

34 

38 

39 

2.491 

2.571 

2.483 

2.568 

1.073 

1.110 

35 

40 

2.648 

2.630 

1.138 



sa 


Table 4.7: Convergence Evaluation of Heuristic 3. 


SI. Problem Ho. of Range of Total 

Ho. Size Problems Ho. of Iterations 

Solved 


Range of Percentage 
Improvements 


11 

o 

6 

r~t 

3 

6 

- 

7 

28.34 


65.34 

2 

>~T 

7 

3 

4 

- 

10 

22.84 

— 

38.85 

3 

8 

4 

4 

- 

8 

9.73 

— 

59.83 

4 

9 

4 

7 

— 

10 

55.27 

— 

83.97 

5 

10 

4 

6 

- 

8 

60.44 

— 

64. 28 

6 

11 

4 

5 

- 

9 

34.00 

— 

88.13 

7 

12 

5 

6 

- 

8 

48.64 


58.13 

3 

13 

4 

6 

— 

9 

77.77 

— 

80.23 

9 

14 

4 

6 

— 

10 

64.66 

— 

89.21 

10 

15 

4 

6 

— 

12 

46.67 

— 

65.15 

11 

16 

3 

7 

- 

9 

38.84 

— 

57.36 

12 

17 

4 

5 

- 

8 

43.86 


68.94 

13 

18 

5 

6 

— 

9 

35.73 

— 

58.75 

14 

19 

3 

8 

- 

9 

28.15 

— 

64.42 

15 

20 

4 

7 

— 

9 

20.15 

— 

52.34 

16 

21 

4 

9 

- 

12 

52.54 

— 

71.67 

17 

22 

3 

7 

- 

10 

32.73 

— 

55.25 

18 

23 

4 

7 

— 

14 

58.63 

— 

84.58 

19 

24 

4 

6 

- 

8 

26.98 

— 

63.72 

20 

25 

4 

6 

— 

14 

27.06 

— 

85.58 

21 

26 

3 

6 

— 

10 

52.93 

— 

88.90 

22 

27 

5 

6 

— 

9 

54.79 

— 

88.63 

23 

28 

4 

6 

— 

8 

29.36 

— 

89.06 

24 

29 

3 

5 

- 

10 

26.63 

— 

86.31 

25 

30 

4 

6 

- 

8 

29 .11 

- 

62.46 

26 

31 

4 

5 

— 

9 

26 .38 

— 

64.76 

27 

32 

4 

5 

- 

8 

9.22 

- 

35.27 

28 

33 

4 

7 

— 

10 

42.82 

- 

64.17 

29 

34 

4 

5 

- 

8 

28.65 

- 

68.48 

30 

35 

3 

6 

— 

9 

33.32 

- 

47.46 

31 

36 

3 

9 

- 

10 

37.40 

- 

53.58 

32 

37 

3 

8 

— 

8 

47.34 

— 

58.93 

33 

38 

4 

9 

- 

9 

50.05 


69.28 

34 

39 

4 

7 

— 

9 

51.04 

- 

73.82 

35 

40 

4 

6 


8 

42.92 

— 

81.34 
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Table 4.8; Range of Percentage Cost Savings Obtained by 
Heuristic One as Compared to the Procedure 
in which All the Items are Ordered at Each 
Set-up. 

(Total Problem = 131) 


SI. Problem Humber of Range of Percentage 

No. Size Problems Cost Savings 

Solved 


1 

6 

5 

52.79 

— 

94.84 

2 

7 

5 

9.88 

— 

85.52 

3 

8 

4 

22.31 

— 

85.27 

4 

9 

4 

57.32 

— 

75.53 

5 

10 

5 

28.72 

— 

94 . 76 

6 

11 

4 

11.72 

— 

72.87 

7 

12 

4 

19.94 

- 

76.93 

8 

13 

5 

38.32 

— 

84.50 

9 

14 

5 

16.99 

• 

87.73 

10 

15 

3 

58.20 

- 

82.47 

11 

16 

4 

67.37 

- 

94.08 

12 

17 

4 

68.62 

— 

85.68 

13 

18 

4 

32.64 

— 

87.50 

14 

19 

5 

77.70 

- 

92.08 

15 

20 

3 

82.33 

- 

85.94 

16 

21 

4 

81.27 

- 

82.82 

17 

22 

4 

48.82 


79.03 

18 

23 

3 

32.74 

- 

44 . 25 

19 

24 

3 

28.35 

- 

60. 26 

20 

25 

5 

40.82 

- 

62.35 

21 

26 

3 

70.76 

— 

97.08 

22 

27 

4 

44.26 

- 

76 . 56 

23 

28 

4 

21.79 

- 

53 . 74 

24 

29 

4 

52.35 

— 

74.86 

25 

30 

3 

21.59 

- 

42.73 

26 

31 

3 

20.95 

- 

68.74 

27 

32 

3 

4.05 

- 

29.87 

28 

33 

3 

43.68 

- 

68.39 

29 

34 

3 

22.15 

- 

24.12 

30 

35 

3 

29.22 

- 

52.19 

31 

36 

4 

18.04 


19.84 

32 

37 

3 

31.48 

- 

33.61 

33 

38 

4 

28.03 

- 

31.17 

34 

39 

3 

34.06 

- 

37.90 

35 

40 

3 

31.96 

- 

64.50 
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